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Abstract The relativistic amplitudes of pion photoproduction are evaluated by dispersion relations at 
t — const. The imaginary parts of the amplitudes are taken from the MAID model covering the absorption 
spectrum up to center-of-mass energies W = 2.2 GeV. For sub-threshold kinematics the amplitudes are 
expanded in powers of the two independent variables v and t related to energy and momentum transfer. 
Subtracting the loop corrections from this power series allows one to determine the counter terms of covari- 
ant baryon chiral perturbation theory. The proposed continuation of the amplitudes into the unphysical 
region provides a unique framework to derive the low-energy constants to any given order as well as an 
estimate of the higher order terms by global properties of the absorption spectrum. 
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^ 1 Introduction 



, In a recent contribution we studied the Fubini-Furlan- 
Rossetti (FFR) sum rule, which relates the anomalous 
^ magnetic moment to single-pion photoproduction on the 
. nucleon pp. This sum rule was derived on the basis of 
' current algebra and PCAC in the chiral limit of mass- 
, less pions 0. It requires a continuation of the production 
' amplitude to the threshold kinematics of massless pions, 
' which is of course outside of the physical region. We eval- 
uated the FFR sum rule by dispersion relations (DRs) 
at t = const., using the imaginary parts of the MAID 
model U as input for the dispersion integrals. As was 
to be expected, we obtained large corrections to the sum 
rule in the physical threshold region due to the finite pion 
mass. However, the FFR sum rule was found to be closely 
obeyed by the production amplitude in the unphysical re- 
gion close to the threshold for massless pions. 

Our further discussion is facilitated by a look at the 
Mandelstam plane shown in fig.^for the isovector photon. 
The two independent kinematical variables are chosen to 
be v and t, related to energy and momentum transfer, 
respectively (see section |21 for a detailed discussion). The 
physical region lies between the solid lines labeled 9 = 
(forward scattering) and 9 = 180° (backward scattering), 
which meet at the production threshold described hy v = 
fc'thr ~ and t = tthr ~ where is the pion 

mass. In the case of a massless pion the threshold moves 
to the origin of the Mandelstam plane, and it is the pion 
production amplitude at this point that is related to the 
anomalous magnetic moment by the FFR sum rule. 



The triangle near the origin is the region where no 
production can take place. It is bounded by the dotted 
lines s = {Mm + Mt^Y and, for the isovector photon, 
t = 9M^. In the former case, the variable s indicates 
the square of the center-of-mass (cm.) energy in the "s- 
channel" + N ^ tt -\- N , which has to be large enough to 
produce a pion and a nucleon with mass M^. The other 
boundary refers to the crossed or "t-channel" 7+7r ^ NN 
whose lowest inelasticity is due to the reaction 7 -I- tt ^ Stt 
requiring a cm. energy of at least three pion masses. 

Whereas the production amplitudes in the described 
triangle can be defined as real functions, they become com- 
plex outside of the triangles due to rescattering or com- 
peting reactions, and therefore the boundaries are lines of 
singularities for these functions. Further singularities ap- 
pear for intermediate one-nucleon and one-pion states in 
the form of pole terms, e.g., s = Af^ and t = AfJ. Since 
these pole terms are well known we subtract them from 
the full amplitudes and refer to the remaining terms as the 
dispersive amplitudes. These amplitudes can then be ex- 
panded in a real power series in v and t, with a convergence 
radius given by the onset of pion-nucleon rescattering or 
the opening of two- and three-pion t-channel reactions for 
the isoscalar and isovector photon, respectively. 

It is the aim of our present paper to derive the 
power series for the pion photoproduction amplitudes with 
MAID05 as an input for the imaginary parts, and thus to 
provide a framework that determines the low-energy con- 
stants (LECs) of chiral effective field theories by global 
properties of the nucleon's excitation spectrum. 
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The general structure of the relativistic ampUtudes was 
already discussed in the early 1990's in the context of a 
then existing puzzle for neutral pion photoproduction at 
threshold @]. By explicitly calculating the pion loop di- 
agrams, Bernard et al. [S] showed that the Taylor coeffi- 
cients of the expansion in v and t were divergent in the 
limit of massless pions, which invalidated the proofs for a 
low-energy theorem that was based on non-singular coefh- 
cients in the chiral limit. Since there appeared to be com- 
plications in the power counting for the relativistic theory, 
the problem was then reformulated in heavy baryon chi- 
ral perturbation theory (HBChPT), which allowed for a 
strict correspondence between the loop expansion and the 
expansion in small external momenta and quark (or pion) 
masses. Being based on a 1/Mm expansion, the frame- 
work of HBChPT leads however to shifts of the pole po- 
sitions, with the result of a spurious behavior of the am- 
plitudes in the unphysical region near the origin of the 
Mandelstam plane pp. These shortcomings of HBChPT 
have been noted often before, and the newly developed 
manifestly Lorentz-invariant renormalization schemes [H] 
now provide a covariant treatment as well as a consistent 
ordering scheme. In particular Bernard et al. T have re- 
cently analyzed the FFR sum rule in the framework of 
infrared regularization of BChPT and obtained both a 
good agreement with the threshold data and the expected 
smooth V dependence of the amplitudes in the (unphysi- 
cal) sub-threshold region. 

We proceed in section |21 by defining the kinematical 
variables and recalling the invariant amplitudes and their 
dependence on the multipoles. In the following section 
we discuss the convergence of DRs at t = const and show 
that the extrapolation of the amplitudes into the region 
near v = t = Q should not be a problem. In section 4 
we discuss the special role of t-channel vector meson ex- 
change and the high-energy contribution to the dispersion 
integral. Our results are then presented in section to- 
gether with a comparison to the data and to covariant 
BChPT. In sectional we close with a short summary and 
an outlook. 



The physical s-channel region is shown in fig.^ Its upper 
and lower boundaries are given by the scattering angles 
9 — and 9 — 180°, respectively. The pion and nucleon 
poles lie in the unphysical region on the straight lines t = 
Ml (pion pole), Vs = (s-channel nucleon pole), and 
Vu = (u-channcl nucleon pole), where 



4Mn 



(4) 



We note that all the involved particles are on their mass 
shell at the point (i' — 0, t — M^). 

The threshold for pion photoproduction lies at 



t'thr 



tthv 



MIMn 



Mn + M^ 

In the pion-nuclcon cm. system, we have 
k^ = {\k\,k), q^^ = {u;,q), 



(5) 



(6) 
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(7) 



2 Kinematics and invariant amplitudes 

Let us first define the kinematics of pion photoproduction 
on a nucleon, the reaction 

j{k) + N{p,)^7r{q)+N'ipf), 

where the variables in brackets denote the four-momenta 
of the participating particles. The familiar Mandelstam 
variables are 

s^{p, + k)l t={q-k)l u={p,-q)l (1) 



and 

v={s-u)/AMn (2) 

is the crossing symmetrical variable. This variable is re- 
lated to the photon lab energy E!:^'' by 

t-Mi 



v = E: 



lab 



AM, 



N 



(3) 



with W — -v/s the cm. energy. 

The transition current operator for pion photoproduc- 
tion can be expressed in terms of 4 invariant amplitudes 

HSl, 

J^=^A.(^,t)Mf, (8) 



with the four- vectors M^ given by 

Mf = -1*75(7^^-^7^) , 
M!^ = 2^75 (Pf^k-q-q^'k-P) , 
M^ = -IJ5 {n^k-q~ Ijtq^') , 

M^ =. -2*75 (7^ fc • F - ^P^) - 2Mn Mf , (9) 

where P'^ — (j)^ +Pj)/2 and the gamma matrices are 
defined as in Ref. |lT)j . 



B. Pasquini et al.: Dispersion Theory and the Low Energy Constants for Neutral Pion Photoproduction 



3 



The invariant amplitudes can be further decomposed 
into three isospin channels Aj {I = +, 0, — ), 

4° = ^|"'ie"3V + + Ai+'<5a3, (10) 

where are the Pauli matrices in isospace. The physical 
photoproduction amplitudes are then obtained from the 
following linear combinations: 











Ai{jn - 









l(0) 



,(0) 



(11) 



The FFR sum rule is derived from the neutral pion photo- 
production amplitude in the limit of g'^ ^ 0. As we note 
from eq. the four- vectors M^, M^*, and M4 vanish, 
and only the four- vector survives in that limit. 

Corresponding to their behavior under crossing, the 

amplitudes ^i^'4'' and ^3""* are even functions of v and 
satisfy a DR of the type 



ReAfiv, t) 



4'P°^°(z.,t) + -7' / du 



,v' lxnAl{v',t) 



(12) 



whereas the amplitudes ^3^''^'' and A\ 2 a ^'^s odd and 
therefore fulfil the relation 

ReAl{v,t) = 



A 



I ,polc 



(-,0 + -^ 



^^,ImA[(^^ (13) 



The nucleon and pion pole contributions are given by 



4 /,po/e 



A I , pole 
^2 



s ~ M% u - Ml 



eg-nN I 1 

2 



3 ~ 2Mn 2 \s~Ml u~Ml 



A 



I, pole _ effTi-jV K 

~ 2Mn 2 



1 



s - Ml ^ u-M} 



, (14) 



— = 1, = Kp — Kn, and 



with = e" — — e 

where Kp and k„ are the anomalous mag- 
netic moments of the proton and the neutron, respectively. 
Additional pole contributions from the t-channel vector 
meson exchange are discussed in section 0] 

The covariant amplitudes Ai can be expressed by the 
CGLN amplitudes J^i jSgOj as follows: 

, M^it - Ml) , Mj,(Ef + Ma,) {t - Ml)_ 



{W-Mj^Y q 



- Ml 



(15) 
(16) 



W + Mn 
W-Mn 



2{W-Mf,)) q 



t-Ml 



2{W + Mn) J 



iEf+MN)^f, 



(17) 



W - Mn [ q 



t - Mi J^. 



t - m; 



2(T^-MAr) q 2(1^ + Mat) 



(18) 



where TV = Att/^{E., + Ma^) {Ef + Mn), q = \q\ and 
all variables are expressed in the cm. frame. Below the 
Z\(1232) resonance, we may limit ourselves to the S- 
wave multipole i?o+ and to the three P-wave multipoles 
Mi+, and Mi_. In this approximation, the CGLN 
amplitudes take the form 



Ti Eo+ + 3(7\fi+ + E^+ ) cos 61 , 
T2 2Mi+ -I- Ml- , 

3iEi+ - Mi+) , .5^4 ^ , 



(19) 



where 9 is the cm. scattering angle, which is related to 
the Mandelstam variables by 

The P-wave contributions are often expressed by the three 
combinations 



Pi = 3£;i+ + Mi+ - Mi_ 
P2 = 3Ei+ - Mi+ -I- Mi_ 
P3 = 2Mi+ + Afi_ . 



(21) 



With these definitions the multipole expansion of eqs. 
(I15|l - (I18f) can be cast into the form 

^,W + Mn f W{t-M'i)\ - 



W - Mn 
MNit - Ml) 



1^2 _ Ml 



P2 



t 



Ao 



W -M. 



N 



W ^Mn 

{P2 - P3 + . . .} , 



- Ml 



P3 + 



(22) 
(23) 
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^3 = 777 77— ^0+ + + -777^ 7T^ I ^1 



W-M, 



N 



A, 



2{W + Mm 



(24) 



0+ 



W{t-Ml) \ [t-Ml) 
W^-M% ) 2{W - Mm) ' 
t - 

W + Mn + 77777^^^ ) ^3 



2{W + Mn) 



(25) 



with Pi — Pi/q and the ellipses denoting the higher partial 
waves. 

The kinematical factors simplify at threshold, and 
eqs. (|22|) - (|25|l take the exact form 

^l(i^thr, ithr) = Mhr 

X <^ Eo+ - — — Pa 



^2(^'thr, ^thr) 



1 + M (l + A*)(2 + /i) 



(Pg + 6A/JV i?) I , 
(26) 



MA,(2 + /i) 



P2-P3- 6MnD 



(27) 



4^ / ^ -^thr f ^ (2 + m)2M^- 

- M^(2 + M) I °+ + 2(1 + ^) 
fiM^ ^ 3(2 + ^)MArAf^ ^ 



2(1 + 



1 + M 



(28) 



- Af^(2 + M) r°+ " 2(1 + 
(4 + 6M + /i^)A% - _ 3M^ -1 
2(1 + Ai) ' l + A^ J' 



(29) 



where = 47rVTT7^/Af^, ^ = M^/Mn, D = (Af2+ - 
E2+ — M2- — E2-)/q^, and all the multipoles have to 
be evaluated at g = 0. In particular we note that the 
amplitude Pi does not appear in eqs. because 
its kinematical prefactor vanishes at threshold. 

Having determined the invariant amplitudes Ai, we 
can combine these results to construct the CGLN ampli- 
tudes J-i as follows: 



= J{Ei + MN){Ef + MN) 



W-M. 



N 



8nW 



f - M 2 

x^Ar + (H^-M.)A.^ ^^^_^;^^ (A,- A,) 



(30) 



/ g, -Mn W + Mn 
^ \lEf + MN SttW ^ 

f — M ^ 

x< - Ai + {W + Mn) Ai - ——^^ {As - A^) \ , 



2{W + Mn) 



(31) 



^3 = ^{E,-MN){Ef+MN) 

x\{W-Mn)A2+A3-A4 



W + Mn 
8t:W 



(32) 



^ / + Mn W-Mn 2 
\lEf + MN 8ttW ^ 

x\ -{W + Mn) A2+A3- Ai 



(33) 



3 Convergence of dispersion relations at 
t=const 

We recall that the integration path for DRs dX t = const 
is fully contained in the physical region only in the special 
case oit — tthr, which path passes through the threshold 
of pion photoproduction (y = i^thr, t = ^thr)- bi all other 
cases the dispersion integrals run through an unphysical 
part of the Mandelstam plane between s = {Mn + Af^)^ 
and the physical threshold value for s, which corresponds 
to forward scattering {6 = 0) for t > tthr and backward 
scattering {9 = 180°) for t < ithr- The only successful 
procedure to construct Im Ai{s,t) in the unphysical re- 
gion has been the multipole expansion of the CGLN am- 
plitudes and the subsequent insertion of hiiJ-'i{s,t) into 
eqs. (|15|I - H18(I . What about the convergence of this expan- 
sion? According to Mandelstam ^] the scattering ampli- 
tudes can be represented by the pole terms and a dou- 
ble integral over the spectral regions A^u-, Ats, and Agt 
shown in fig. ^ The boundaries of these regions result 
from the lowest mass intermediate states possible in any 
pair of the variables s, t, and u. We hasten to add that 
"possible" refers to solutions with real values of s, t, and 
u and all particles on their mass shell, however there is 
no overlap between the double spectral regions and the 
physical regions. Therefore the Mandelstam representa- 
tion has never had any direct consequences for the data 
analysis. This representation has however important con- 
sequences because it reflects maximum analyticity, in the 
sense that the only singularities are given by the poles 
due to one-particle intermediate states and the cuts due 
to the onset of particle production channels. In particu- 
lar, one-dimensional DRs such as DRs at i = const can 
be straightforwardly derived from the Mandelstam repre- 
sentation 12 . After subtraction of the pole terms from 
the full amplitudes Ai , the range of convergence is deter- 
mined by the nearest singularity in the Mandelstam plane, 
which arrives when the line t = const becomes tangent to 
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the double spectral region In the case of the multi- 
pole expansion, the convergence is based on the following 
mathematical lemma: If a function /(z = x + iy) is ana- 
lytic inside and on an ellipse C whose foci are at the points 
{x = cos 9 = ±1, y = 0), it can be expanded in a Legendre 
series for all points in the interior of the ellipse C |14j . 

The relevant relations for t may be obtained from 
eq. (|20|l and the kinematics of eq. (TJ, 



The first of these curves defines the upper limit of conver- 
gence by the asymptote iuppci(s — ^ oo) = 4M^. The lower 
limit is again obtained by reflection according to eq. H38(l . 
and its maximum at s sa 1.65 GeV'^ yields the lower limit 
of convergence at t„iin = —0.525 GeV^. Altogether then 
we find that DRs should converge in the following strip of 
the v ~ t plane: 



t{s) = Ml- 2w(s)fc(s) + 2q{s)k{s) cos ( 



(34) 



The center of the ellipse lies on the line 9 = 90°, 

to{s) ^ Ml - 2oj{s)k{s) , (35) 
and the foci correspond to forward and backward scatter- 



ing, 



t±{s) =to{s)±2q{s)k{s). 



(36) 



The ellipse of convergence can now be stretched until its 
upper or lower tangent t — const touches the nearest dou- 
ble spectral region of the Mandelstam representation. For 
the isovector amplitudes the boundary of the double spec- 
tral region a[^^ is given by 



iupper(s) = dMl + 



8M4 (3s - Ml, + Ml) 



s - (Mjv + Af,)2] [s - {Mn - M,)2] ' 

(37) 



with the asymptote tuppcr(s — oo) = QMl as the upper 
limit of convergence of DRs at t = const. The lower limit 
is given by the reflection of eq. I|37|) at the center of the 
ellipse, eq. H35|) . This leads to 

tlowcr(s) = 2io(s) - iuppor(s) ■ (38) 

The maximum of this function lies at s ~ 1.32 GeV^ and 
takes the values i^in = -0.376 GeV^ and -0.352 GeV^ 
for charged and neutral pions, respectively. 

The corresponding boundary for the isoscalar ampli- 
tudes j4^(' is defined by the lower value of the two curves 



tl:lpeM^Q^{s) + ^Qns)+R,{s), t=l 



and 2 , 



(39) 



with 

(.) = 2M' + (9. + 31M2, - 28Ml)Mt 

' " - {Mm + 2Af,)2] [s - {Mm - 2M,)^] ' 

(40) 

O (,) = 8M2 + ^{9s + Ml,-Ml)M^ 

' ^ [s- {Mm + M.)2][s - {Mm - M.)2] ' 

(41) 

J. , . (4M2 - Ml)Ml 

^'^'^ \s - {Mm + 2M,)2][s - {Mm - 2M^Y\ ' 

(42) 

^ , ^ _ 16(M2, - Ml)Ml 

' Is- {Mm + M.)2] [s - {Mm - M.)2] ' 

(43) 



Af"^ : -0.376 GeV^ < t < 0.175 GeV^ 
: -0.525 GeV^ < t < 0.073 GeV^ , 



(44) 



For completeness we recall two other critical values below 
these limits. With decreasing t the s- and u-channel cuts 
approach and touch at i/ = 0, i = — 2M^(2MAr + M^r) « 
—0.564 GeV^. Up to this point, and with some care also 
in the region below, DRs should still be valid, except that 
the Legendre expansion is no longer convergent. The final 
break-down of DRs occurs if t — const is tangent to the 
double spectral region Ait'^'' given by the relation 

[s _ (Afjv + M^)^] [s - {Mm - M^)^] [u - {Mm + M^)^] 

x[u- {Mm - M,)2] - (4A/4 - MI)MI 

X [2su - 2(Af^ - Ml){s + u) + 2M% ^ M^] = . (45) 



This happens at = and t = —1.10 GeV^ and 
— 1.06 GeV^ for charged and neutral pions, respectively. 

For completeness we also mention the work of Oehme 
and Taylor |12j . On the basis of relativistic quantum field 
theory, especially causality and the spectral conditions, 
these authors have rigorously proved that a polynomial 
expansion in cos9 converges for —12Ml < t < 0, which 
yields a considerably smaller range than in the case of the 
Legendre expansion of Ref . jEl • 



4 Vector mesons and Regge tails 

The t-channel exchange of vector mesons plays an im- 
portant role in neutral pion photoproduction. The vector 
meson contributions to the invariant amplitudes take the 
form: 



AY{t) = 



t 



2MmM^ t-mlr' 



A^{t) = - 



2MnM^ t- 



(T) 



A^{t) 



0, 



e\v gy 
M^ 



(46) 



where Ay denotes the coupling of the vector meson {V ~ 

uj,p) to the 77r'^ system and g\)^''^^ its vector or tensor 
coupling to the nucleon. Due to their isospin structures, 
the uj contributes to the isospin amplitude I = + and the 
p to / = 0. In the 2005 version of MAID the coupling 
constants for the w (782) (in bracket: p(770)) are: Ay = 
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0.314(0.103), g^y^ = 16.3(1.78), g^P = -15.4(22.6). 
Whereas the coupUngs Ay for the V "fir process are 
essentially known, the vector and tensor couplings gy ' 
to the nucleon are less well determined. However, all anal- 
yses agree that the lu exchange is quite essential for neutral 
pion photoproduction, whereas the p exchange is almost 
negligible. 

In the zero-width approximation the vector meson 
poles at t = rriy should play a similar role as the nu- 
cleon poles at = ±j^B(t) and, in the case of charged 
pion production, the pion pole at t = . All these pole 
terms can not be obtained by the dispersion integrals of 
eqs. (fT^ and (fT!I|l but have to be added to the dispersive 
contributions. Of course, one has to make sure that the 
pole terms are not modified by errors of numerical origin 
and by the derivation of the absorptive amplitudes from 
the data. In order to eliminate such double-counting in the 
case of charged pion electroproduction, von Gehlen 
has proposed to subtract the term 



physical region by means of Legendre polynomials breaks 
down for t < —19 AI^ (see section OJ. 

In particular we note that the t-channel exchange of 
vector mesons does not contribute to the crossing-odd am- 
plitudes A''^'^\ which can only receive pole contributions 
from the exchange of axial vector mesons. 

At energies above the resonance region, the Regge 
model provides a convenient description of the cross sec- 
tions. It it obtained by replacing the p and oj propagators 



1 



v' lim. 



[f ■ 



Af2 )ImAf(^^',0 



(47) 

from the dispersion integral of eq. I|12|l for the appropri- 
ate electroproduction amplitude (i = 5, / = — ). In the 
same spirit we have checked whether the dispersion inte- 
gral can provide a pole structure at i = rn^. The result is 
negative. Even though the vector meson background plays 
an important role in the unitarization process of MAID, 
there is no indication of a vector meson pole term similar 
to eq. H47|l in our calculation. It is even more surprising 
that the dispersion integrals for the threshold amplitudes 
change only by a few percent if we drop the vector mesons 
in the construction of the absorptive amplitude, whereas 
the vector mesons yield 20 % and 50 % of the dispersive 
contributions for A2 and A4, respectively. 

In view of the problem to reproduce the vector me- 
son poles by the dispersion integrals of eq. (|12l) , we accept 
the Mandelstam hypothesis JTj that the amplitudes are 
the sum of all pole terms plus two-dimensional integrals 
over the double spectral region. The one-dimensional DRs, 
e.g. aX t = const, follow from this representation, as has 
been proved for pion photoproduction by Ball JJ5,. Al- 
ternatively, we could subtract the DRs at = 0, which 
introduces an unknown function Al (0, t). This function is 
real in the region of small t, and in principle can be con- 
structed from its imaginary part |17) by an integral along 
the t-axis (see fig. DJ. For i > the absorptive amplitude 
is dominated by the reactions 7 -I- tt ^ 27r (/ = 0) and 
Stt (J = 4-) whose resonant parts are dominated by vector 
meson production occurring aX t m rriy « 34M^. How- 
ever, if we follow the negative t-axis to t = —29M^ and 
below we also find absorptive amplitudes, starting with 
S-wave production, followed by P-wave contribution from 
Z\(1232) excitation, and finally effects from the totally 
unknown double spectral region for t < —58M^. Unfor- 
tunately, these contributions at t < can not be derived 
from the data basis, because the extrapolation to the un- 



{t — rriy) ^ by the Regge propagators |18II19| 

with av{t) = aQ+ait describing the Regge trajectory and 
So a free parameter. Typical values for the w trajectory 
(in brackets: p trajectory) are ao = 0.44 (0.55) and ai = 
0.9 GeV-2 (0.8 GeV-^). As a result the Regge amplitudes 
have a typical l/y/s = 1/W behavior for small values of 
t. 

We have estimated the high-energy tails of the dis- 
persion integrals by use of the following models for the 
imaginary parts of the amplitudes: 

(I) a 1/W^ tail fitted to MAID in the energy range of 1.8 - 
2.5 GeV, 

(II) an energy dependence according to eq. fitted as 
above, and 

(HI) phenomenological parameterizations of Regge trajec- 
tories and cuts fitted directly to older high-energy 
data pnil?T| . 

While none of these descriptions is completely satisfactory, 
such studies provide a reasonable estimate of the high- 
energy contributions. 



5 Results 

The dispersive contributions to the real part of the ampli- 
tudes Ai{i',t) are determined by the integrals of eqs. (|12|l 
and H13|l . The following figures show the respective inte- 
grands for t = ithr and the two values — and — vthi-. 
The dispersive amplitudes are then obtained by integra- 
tion and multiplication by the factors 2/7r or 2i^/7r in 
front of the integrals. At this point we remind the reader 
that the amplitudes Ai have different dimensions due to 
the traditional definitions of eqs. |(HJ and 0, and as a 
result the integrands are given in different units. We fi- 
nally note that all the following calculations are obtained 
with MAID05, which now describes the amplitudes up to 
W = 2.2 GeV. 

MAID is a unitary isobar model for photo- and elec- 
troproduction of pions in the resonance region. It is con- 
structed from a field-theoretical background with nucleon 
Born terms and t-channel vector meson exchange pole 
terms, both unitarized in K-matrix approximation pil22| . 
The resonance sector is modeled with s-channel nucleon 
resonance excitations in Breit-Wigner form for all four- 
star resonances below W — 2 GeV. The hadronic p and oj 
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coupling constants of the background and the electromag- 
netic resonance couplings ^3/27 S'1/2 that describe 
electric, magnetic and Coulomb multipoles are fitted to 
the world data base of pion photo- and electroproduction. 
Hadronic resonance parameters are taken from the Parti- 
cle Data Tables Latest versions of MAID (MAID03 
and MAID05) describe very well the data in the kine- 
matical range of Wthr < W < 2 GeV and < < 
5 GeV^ |24II25| . Over a wide energy range up to the sec- 
ond resonance region MAID is generally consistent with 
dispersion relations at t = const |26j . 

Figure |2] shows the integrands for the four isospin (+) 
amplitudes. The dashed lines give the integrands a,t v — 0, 
which are strongly dominated by the Z\(1232) resonance 
with a much smaller contribution of the second resonance 
region for Ai and A2. In addition, there appears a strong 
background of S-wave pion production in the case of yli . If 
v approaches the onset of S-wave production (see the full 
line!), the integrand increases dramatically and in prin- 
ciple runs into an inverse square root singularity, which 
eventually leads to the cusp effect clearly seen in tt" pro- 
duction at the tt"*" threshold. The cusp effect is not present 
for A2, because according to eq. (I22|l the S wave does not 
contribute to this amplitude, and S waves are also strongly 
suppressed in the case of A3 and A4 due to kinematical 
enhancement of the P-wave combinations P2 and P3, re- 
spectively. 

The integrands for the amplitudes A^^'' are displayed 
in the following fig. O Since the isoscalar photon can not 
excite the A resonance, we now find a competition of con- 
tributions from the second and third resonance regions 
plus a large cusp effect for A^^\ Altogether these inte- 
grands are considerably smaller than in the case of the 

The presented integrands show a reasonable conver- 
gence for the higher energies. In fig. 01 we investigate the 
question of convergence further by decomposing the imag- 
inary parts of the amplitudes into a partial wave series. 
The figure shows the dispersive part of the n^p am- 

plitudes as function of v and a,t t — ithr- While both S and 
P waves in the imaginary part yield large contributions to 
Al, the other three amplitudes are essentially determined 
by the P waves, with some higher partial wave contribu- 
tions in the case of A2 and A3. 

Since we want to extrapolate the amplitudes into the 
unphysical region at small t values, we also study the dis- 
persive parts of the isospin amplitudes as function of v for 
a series of t values in the range > t > — lOAf^. As 
shown in fig. |S1 the t dependence develops quite a regular 
pattern, with the cusp moving to larger u values with in- 
creasing values of t. At the same time the physical thresh- 
old (see the line 6 — 180° in fig.CJ moves to smaller values 
of u down to the minimum ai v — vthi-, from whereon it 
increases again (see the line 6* = in fig. 

The comparison with the experiment :.27^ in fig. 
shows that the dispersion integral by itself misses the 
threshold data for the amplitudes A2 — A4. If we add the 
t-channel p and ui poles according to MAID05, we obtain 
an almost perfect agreement for Ai, A2 and A4. The ap- 



parent discrepancy between theory and experiment for ^3 
is an open question, which will be discussed later. 

After these tests we feel safe to expand the amplitudes 
about the unphysical point at = and t = where 
all the involved particles are on their mass shell. For this 
purpose we proceed as in our previous work ^ by casting 
the first amplitude in the form: 



AfP(j.,t) = Ai(z.,i)-Ar"^(j.,t) 



(49) 



where k is the anomalous magnetic moment in the re- 
spective isospin or physical channel and Ai{u,t) the di- 
mensionless "FFR discrepancy" . For convenience we also 
express the three other amplitudes by dimensionless func- 
tions Aiiv, t), 



Ar'{'^,t) = ^A2{.,t), 



2M^ 
2M^ 



(50) 



The functions Ai (v, t) are regular near the origin of 
the Mandelstam plane and can be expanded in a power 
series in v and t or vb- As is evident from eq. ;/ is 
0{Mt^) and vb is 0{M^/Mn) in the region of interest. 

Therefore the crossing-even amplitudes ^1^'"' and 
have the expansion 

A{u, t) - 5oo + ^20 y^Ml + <5o2 i^b/M^ + S40 ^VK 

+S22 v^vbImI + <5o4 vllMl + . . . , (51) 



with the lowest expansion parameters given by 

^20 



Z\(0, Ml) , 



2 9j/2 



i/=0 ■ 



d 



002 



4A%A4^Z\(0,t)|^^ 



(52) 



(+.0) 



In the case of the crossing-odd amplitudes Z\; 
a\ 2 4 the corresponding expansion takes the form 



and 



A{u., t) = Sio ly/A'U+Sso ly^M^+Su i^yB/Ml+. . . , (53) 



with the lowest expansion parameter 



<5io = 



d 



(54) 



In order to study the convergence of the multipole series 
involved in evaluating the dispersion integrals of eqs. (|12|l 
and H13|) from the imaginary parts of the amplitudes, 
we compare the results obtained dX t = with val- 
ues obtained at different t. Tabled shows our results for 
A'f'^"\0,t) for t = Ml, 0, and ithr- The total value at 
t = Ml corresponds to 6qq — —0.078. The slow decrease 
with decreasing values of t is due to the term So2i^B- 
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t 


Ml 







L = 


0.571 


0.579 


0.586 


L = 1 


1.209 


1.150 


1.098 


L = 2 


-0.109 


-0.089 


-0.071 


L = 3 


0.044 


0.054 


0.062 


sum 


1.715 


1.695 


1.676 



Table 1. Contributions of the dispersion integral to Hp + 

A[^'"°''{0,t) at t ^ Ml, 0, and tthr for the pn" channel. The 
rows i = to L = 3 indicate the contribution of the S, P, D, 
and F waves to the dispersion integral. 



t 


Ml 





ithr 


L = 


0.11 


0.12 


0.12 


L = 1 


0.21 


0.21 


0.22 


L = 2 


-0.003 


-0.002 


0.000 


L = 3 


-0.000 


0.000 


0.000 


sum 


0.32 


0.33 


0.34 



Table 2. Contributions of the dispersion integral to 

■^^-j-Zij'"' \y,t) for the pvr" channel at fixed v = Q. See 
table □ for further notation. 



t 


Ml 





ithr 


L = 


-0.22 


-0.23 


-0.24 


L = 1 


1.62 


1.66 


1.71 


L = 2 


-0.54 


-0.57 


-0.59 


L = 3 


-0.29 


-0.26 


-0.24 


sum 


0.57 


0.59 


0.63 



Table 3. Contributions of the dispersion integral to 
AMn M^-^A'f"'°\Q,t) for the p7r° channel at fixed u ^ Q. 
See table rrl for further notation. 







(520 


So2 




-0.08 (+0.04) 


0.32 


0.57 (+1.11) 


A2 


-4.54 (-1.87) 


-1.26 


3.47 (-1.55) 




-2.23 






Ai 


13.02 (+8.17) 


2.23 


-2.48 (+6.98) 



Table 4. The leading expansion coefficients for the pn ampli- 
tudes from the dispersion integral (see eqs. H49|l -I)54| | for defini- 
tions) and the vector meson t-channel contributions (in brack- 
ets). 



The following tables |21 and 13 show the t dependence 
of the derivatives with regard to i/"^ and t at fixed i/ = 0. 
The total values a.tt = Ml yield the constants ^^'^ ^02 > 
respectively. The small changes with decreasing values of 
t are due to the higher expansion coefficients S22 and S^^. 
We note that the curvature in v (table is essentially 
determined by the S and P waves which add with a 1:2 
ratio. However, the slope in t (tabic |2J) has considerable 
contributions from the higher partial waves. Although the 
P wave yields by far the largest contribution, the other 
partial waves conspire to reduce the total slope to less 
than 30% of the P-wave result. 



In table 01 we list the leading expansion coefficients of 
the 4 invariant amplitudes for the pn^ channel. The ad- 
ditional contributions of the vector meson poles to the 
coefficients ^oo and 6q2 are given in brackets. It is obvi- 
ous that these i-channel effects play an important role 
in neutral pion photoproduction. The strong competition 
between s- and t-channel contributions to Sqq and ^02 re- 
flects the previous discussion on subtracted DRs. Accord- 
ing to eqs. (|49|I - H51() the subtraction functions Af {v — 0, t) 
are determined by the respective expansion coefficients 
^00 + ^021^3 /My^ + ■ ■ and on the other hand the DRs 
would sample information at v = from both t-channel 
reactions (for t > 0) and s-channel reactions extrapolated 
into the unphysical region (for t < 0). 

The numbers in table 01 should be compared to an 
expansion of the loop plus counter terms in covariant 
BChPT. Such a calculation has been performed in Ref. [7] 
by evaluating the third-order loop corrections and sup- 
plementing them by a fourth-order polynomial. Since the 
fourth-order loop corrections are large, the resulting power 
series is only indicative of the expected LECs j2|. However, 
the coefficients compare favorably with the LECs obtained 
from the earlier HBChPT calculations "SF. Including for 
consistency an additional factor e on the RHS, we obtain 
from eq. (|12() of Ref. |7j the following coefficients from the 
fourth-order polynomial contribution: 6qq = 0, J^o = 0.53, 
= 3.40, S^o = -6.33, Sf^ = -2.58, and 6^^ = 22.4. 
These numbers are in qualitative agreement with our re- 
sults in tabled The differences are due to 

(I) the power series expansion of the loop corrections, 

which has to be added to the polynomial, 
(II) the effects of higher partial waves, particularly with 

regard to the t-dependence given by 602 , and 
(HI) possible s-channel resonances above 2.2 GeV and t- 
channel exchange of heavier objects, not included in 
the dispersive approach. 

As we have shown, the dispersion calculation including 
the vector meson poles is able to reproduce the experimen- 
tal threshold data except for the crossing-odd amplitude 
A3. In order to pin down the origin of this discrepancy, 
we have checked the integrands of the dispersion integrals 
by repeating our calculation with the SAID [55] analysis. 
Whereas the real parts of the MAID and SAID multipoles 
differ in some cases, particularly at the higher energies, the 
imaginary parts generally agree well for W < 2.2 GeV. 
As shown in fig. the dispersive contributions of the 
two models turn out to be quite similar. Whereas the 
experimental threshold values for A2 and A4 can be de- 
scribed after adding the vector meson pole contributions, 
no such remedy exists for the crossing-odd amplitude A3 
for which both SAID and MAID are off by a factor of 
2. We have also checked the high-energy contributions 
by varying the onset of the asymptotic tail in the range 
1.8 GeV^ W < 2.5 GeV and its shape from a simple 1/W 
dependence to various Regge prescriptions. In this way we 
can modify the threshold amphtudes A\'^^, ^l'"', and A*J^^ 
by at most 10 % in the n'^p channel. However, the asymp- 
totic contribution to A*^^ reaches at most 1 % because of 
the better convergence for this crossing-odd amplitude. 
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We have also directly estimated the high-energy tail 
with the parameters given by the Regge models of the 
1970's ^0,21 , which were constructed to fit the data in 
the 5-20 GeV region. Since these models include the ex- 
change of axial vector mesons and Regge cuts correspond- 
ing to many-particle exchange, they also contribute to 
the crossing-odd amplitude ^3. However, the predicted 
strength of the high-energy contribution to A3"' turns 
out to be even smaller and of the order 10^"^ only. On 
the phenomenological level, possible candidates for axial 
vector meson exchange (J^*^ = 1'''"'') are the ai (1260) 
with quantum numbers J*^ = 1^ and the /i (1285) with 
J*^ = 0+. The a\ has the same isospin and G parity as the 
pion, and therefore contributes only for charged pion pro- 
duction. The /i, on the other hand, has positive G parity 
and a branching ratio of about 5 % for the 7p° channel. 
It is therefore a good candidate for p meson photoproduc- 
tion and, together with the ai, also for the Regge tail of 
the helicity-dependent inclusive cross sections as measured 
by the Gerasimov-Drell-Hearn experiment 30^. However, 
there appear to be no clear candidates to solve the 
puzzle by a t-channel pole term. 

Let us finally discuss the multipole content of the rela- 
tivistic amplitudes and the related error bars in the thresh- 
old region. As in our previous analysis we subtract the 
nucleon pole terms, which may vary within about 2 % 
depending on the value chosen for the pion-nucleon cou- 
pling constant g-j^N ■ Since the pole term constitutes about 
85 % of the total threshold amplitudes A|'"' and 
the choice of g^Ar leads to an error of about 12 % in the 
remaining dispersive amplitudes. For Af^'^ and Al'"", how- 
ever, the pole contributions are small and therefore the 
model error of Qt^n can be neglected with regard to the 
dispersive amplitudes. In table[51the threshold amplitudes 
are constructed from the experimental values of the S- and 
P-wave multipoles and the MAID05 value for the D waves. 
The error given in the table is obtained by adding the ex- 
perimental errors for the threshold multipoles. We recall 
at this point that here and in the following all values refer 
to the dispersive contributions only. 

As is evident from table [51 A^'*'' is dominated by the 
S-wave, because the magnetic contributions cancel nearly 
completely. The large S-wave contribution originates from 
the FFR current U and rescattering corrections, which 
are of course included in BChPT by the chirally invariant 
pion-nucleon coupling and the pion loops leading to the 
pronounced cusp effect at the nTr"*" threshold. These ef- 
fects are nicely reproduced by the dispersion integral, and 
the small t-channel pole contribution is well within the ex- 
perimental error bars. More precisely, the vector mesons 
produce large effects of about equal size for both mag- 
netic multipoles, but the discussed cancellation between 
Mi_|_ and Mi__ leads to a reduction by a large factor. 

Quite different physical information is sampled by 
A^^ . Due to the structure of the associated four- vector 
M2, the spin J = 1/2 multipoles i?o+ and Mi_ do not 
appear in this amplitude. It is dominated by Afi+ but 
also receives surprisingly large contribution from the D 
waves whose threshold values are not very well known. In 





Eo+ 




Mi+ 




D 


total 




4.75 


0.08 


2.02 


-1.94 


0.08 


5.00 ±0.25 


A2 





-0.34 


-29.44 





5.13 


-24.7 ±5.9 


A3 


2.36 


-0.65 


-18.67 


14.34 


-0.65 


-3.3 ±1.7 


Ai 


2.36 


0.04 


40.55 


14.34 


0.04 


57.3 ±1.9 



Table 5. The multipole decomposition of the dispersive part 
of the experimental threshold amplitudes, constructed from 
the data of Ref. :27, and the MAID05 value for the D-state 
contributions. The errors include statistical and systematic er- 
rors from the experiment and model errors due to uncertain- 
ties from the pion-nucleon coupling and from the D waves, all 
added in quadrature. 



HBChPT 28; about 90 % of this ampUtude are obtained 
from counter terms. 

In the case of A3'"' we again find a large cancellation 
of the magnetic multipole contributions, which is required 
because the p and lu pole terms have to cancel exactly for 
symmetry reasons. However, the dispersive contributions 
of Mi+ and Mi_ are now large compared to Eo+, and this 
cancellation leads to a large error bar. The amplitude A3 
shows a weak cusp effect, in accordance with the result of 
HBChPT EH]. 

Finally, the amplitude A4'"' is determined by the 
Z\(1232) multipole Mi^ and a remarkably strong Roper 
multipole Mi_. Since no cancellation appears among the 
dispersive contributions, the error bar of this amplitude 
is small. Similarly as in the case of A^^ , the amplitude 
A4'"' is almost totally determined by counter terms in both 
HBChPT ,22, and covariant BChPT (2j. 



6 Summary and Outlook 

We have studied the photoproduction of pions by means 
of dispersion relations at t = const for the relativistic am- 
plitudes Al to A4, using the imaginary parts of the ampli- 
tudes as input to evaluate the real parts. The calculations 
are performed with and compared to the most recent ver- 
sion MAID05, and it is our final goal to put this analysis 
on the basis of dispersion theory. Along these lines the 
present exploratory work on neutral pion photoproduc- 
tion has the aim to (I) test the numerical procedure and 
the data basis by trying to reproduce the precision data 
near threshold and (II) predict the low-energy constants 
of baryon chiral perturbation theory (BChPT) by global 
properties of the excitation spectrum. 

After subtraction of the nucleon and pion pole terms, 
the remaining "dispersive" amplitudes are regular func- 
tions of the Mandelstam variables in the subthreshold re- 
gion. In particular they may be expanded in a power series 
in the two independent variables v and t about the point 
v = 0, t = AI^. This scries converges in a circle with a ra- 
dius determined by the threshold of pion production. The 
singularity at threshold is, of course, well described by the 
loop corrections of BChPT, which also can be expanded 
in a power series below threshold. The difference of the 
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expansions of the full amplitude and the loop contribu- 
tion is described by the low-energy constants necessary to 
apply BChPT to the data analysis. 

We find that the threshold amplitude Al'^^ is well de- 
scribed by the dispersion integral within the experimental 
error. As in our previous work we also obtain a good agree- 
ment with the FFR sum rule that connects the amplitude 
Al to the anomalous magnetic moment of the nucleon. In 
the case of and Af^"^, however, we have to include ad- 
ditional contributions of t-channel vector meson exchange. 
The importance of such effects for neutral pion photopro- 
duction has been noted long ago. As was to be expected, 
our calculations also show that these fixed poles at t = my 
can not be obtained from the dispersion integrals. More- 
over, though these mesons play an important role for the 
unitarization process at the higher energies, their global 
eflFect on the dispersion integrals for the threshold ampli- 
tudes is surprisingly small. If we add the vector meson 
pole terms to the dispersion integrals, we obtain a consid- 
erable improvement also for v4|'"' and A*^^ , and a perfect 
fit would require only a modest change of the coupling 
constants with respect to the MAID values. 

However, the predicted value for Af^^ is at variance 
with the data by more than 2 standard deviations. This is 
a serious problem for the crossing-odd amplitude A3, be- 
cause a t-channel pole contribution requires the exchange 
of an axial vector meson. Unfortunately, the known axial 
vector mesons have either the wrong isospin or the wrong 
G parity for neutral pion photoproduction. 

In view of the apparent discrepancy for Ag'*'" we have 
carefully checked all the ingredients of our calculation. We 
have repeated the integrations using the absorptive am- 
plitudes given by SAID 29 . In spite of occasional differ- 
ences between the MAID and SAID multipole analyses, 
particularly at energies W >2 GeV, both models predict 
quite similar amplitudes up to the Z\(1232) resonance. We 
have further studied the influence of the high-energy re- 
gion by varying the upper limit of integration in the range 
of 1.8 GeV< W ^ 2.5 GeV and by fitting the high-energy 
tail with various shapes from a simple l/W behavior to 
several Regge forms. In this way the threshold amplitudes 
change by about 5 % for A*'"' and a few per cent for ^2^'' 
and A*^^ , whereas the possible error for A3'"' is less than 
1 %. If we directly apply various Regge models, including 
phenomenological Regge poles and cuts, the high-energy 
contribution turns out to be even smaller, particularly in 
the case of A3'"' . Even independent of any model assump- 
tion, the crossing-odd amplitude A3 has the best conver- 
gence property of all the amplitudes, and therefore an ex- 
tremely strong asymptotic contribution would be neces- 
sary to change A3'"' by a factor of 2. 

As we have shown, the experimental error for A|'*'' is 
very large due to a near total cancellation between the 
leading magnetic multipoles. For this reason the D state 
contribution may give rise to some concern. Its contribu- 
tion has been estimated to account for about 20 % of A3'"' , 
which assumption is based on a reasonable extrapolation 
to the threshold region. However, local fits in certain en- 
ergy bins often yield large fluctuations of the D-state back- 



ground. Therefore a large model error for this background 
can not be excluded and may be partially responsible for 
the A3 problem. 

With all these caveats in mind we have evaluated the 
dispersive and vector meson contributions of the relativis- 
tic amplitudes and expanded the result in a power series 
in v and t as discussed in detail in the previous section. 
The results nicely compare with the low-energy constants 
derived from a fit of covariant BChPT to the threshold 
data. Except for the A3 problem, the differences between 
the two approaches are due to additional loop terms. In 
view of the unexpectedly large loop corrections at 4*** or- 
der HBChPT, an extension of the existing 3'''* order co- 
variant BChPT to the next order will be of great general 
interest. 

On the side of the dispersive approach, the simple ad- 
dition of the vector meson pole terms to the dispersive am- 
plitude is of course only justified near threshold where the 
scattering phases are small. In general it will be necessary 
to unitarize the full amplitude comprising the nucleon, 
pion, and vector meson pole terms as well as the com- 
plex contribution of the continuum. This will require an 
iterative procedure including modifications of the model 
parameters encoded in the hclicity amplitudes of the reso- 
nances. A further challenge is to describe the gradual tran- 
sition from vector meson poles near threshold to Regge 
propagators at large energies, without introducing spuri- 
ous singularities and respecting the crossing symmetry of 
the invariant amplitudes. 

In view of the large D-wave background and the deli- 
cate cancellation of various multipoles in A3'"', dedicated 
experiments to measure the D-wave contribution in the re- 
gion between threshold and the Z\(1232) resonance would 
be extremely helpful. Whereas the difference of recoil and 
target polarization in forward direction is directly pro- 
portional to the amplitude A3 in the high-energy limit, 
this observable does not appear to be very sensitive to A3 
in the threshold region. We are presently studying sev- 
eral double-polarization observables with regard to an en- 
hanced sensitivity for A3 and the D-wave background. 
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Figure 1. The Mandelstam plane for the isospin / = + as 
function of f and t. Solid lines: the boundaries of the physical 
region between forward {9 = 0) and backward (6 = 180°) scat- 
tering. Dashed lines: the boundaries of the st and su spectral 
regions. The dashed-dotted line is obtained by reflecting the st 
boundary at the dotted line 8 = 90°. The other dotted lines in- 
dicate the asymptotic boundaries of the spectral regions. The 
left half-plane, which is not shown here, is determined by the 
crossing symmetry, v or s <-> See text for further 

explanation. 
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Figure 2. The integrands of eqs. 11211 and 11311 for the am- 
plitudes A^^' obtained with u = fthr (solid lines) and i/ = 
(dashed lines). 
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Figure 3. The integrands of eqs. 1121 and 1131 for the am- 
phtudes A^"^ obtained with ly — fthr (solid lines) and u = 
(dashed lines). 
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Figure 4. The real parts of the amplitudes A^^^ ' for the 
reaction 7p ^ vr^p as function of v and at t = tthi- Solid lines: 
Dispersive contributions according to eqs. (I12II and 11311 as ob- 
tained with the imaginary amplitudes from MAID05 contain- 
ing partial waves up to Lmax ~ 3. Short-dashed lines: S-wave 
contributions, dotted lines: P-wave contributions, long-dashed 
lines: Sum of D- and F-wave contributions. The dashed-dotted 
lines starting at pion threshold show the real parts taken di- 
rectly from MAID05. 
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Figure 5. The real parts of the amplitudes A^^^ as function of 
1/. The dispersive contributions according to eqs. 1121 and 11311 

as evaluated with the imaginary amplitudes from MAID05 for 
the following values of t: (solid lines), tthr (dashed lines), 
-4M^ (dashed-dotted lines), and -lOM^ (dotted lines). 
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Figure 6. The real parts of the amplitudes A^^^ for the 
reaction ■yp 7r°p as function of u and at t = tthr- Solid lines: 
dispersive contributions according to eqs. I12II and as eval- 
uated with the imaginary amplitudes from MAID05 containing 
partial waves up to Lmax ~ 3. Dashed lines; same results cal- 
culated with SAID. The dashed-dotted lines are obtained by 
adding the vector-meson contributions of eq. 14611 to the MAID 
result. The data points near threshold are derived from the ex- 
perimental values of Ref. |27| for the S and P waves plus the 
MAID correction for the D waves. 



